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ABSTRACT 

Recent developments in numerical relativity have made it possible to follow 
reliably the coalescence of two black holes from near the innermost stable circular 
orbit to final ringdown. This opens up a wide variety of exciting astrophysical 
applications of these simulations. Chief among these is the net kick received 
when two unequal mass or spinning black holes merge. The magnitude of this 
kick has bearing on the production and growth of supermassive black holes during 
the epoch of structure formation, and on the retention of black holes in stellar 
clusters. Here we report the first accurate numerical calculation of this kick, 
for two nonsp inning black holes in a 1.5:1 mass ratio, which is expected based 
on analytic considerations to give a significant fraction of the maximum possible 
recoil. We have performed multiple runs with different initial separations, orbital 
angular momenta, resolutions, extraction radii, and gauges. The full range of 
our kick speeds is 86-116 km s _1 , and the most reliable runs give kicks between 
86 and 97 km s -1 . This is intermediate between the estimates from two recent 
post-Newtonian analyses and suggests that at redshifts z > 10, halos with masses 
< 10 9 M Q will have difficulty retaining coalesced black holes after major mergers. 

Subject headings: black hole physics — gravitational waves — relativity — cos- 
mology: theory 



1. Introduction 

When two black holes merge, the gravitational waves they produce will carry away 
net linear momentum, barring conditions of special symmetry (e.g., two equal- mass non- 
spinning black holes receive no kick). The magnitude of the resulting recoil is important 
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However, it has been shown that almost all of the recoil occurs in the strong gravity regime, 
inside the innermost stable circular orbit (ISCO). This is precisely where analytical treat- 
ments are least reliable. An accurate estimate of the recoil kick therefore requires full numer- 
ical simulations of the final phase of the inspiral, merger, and ringdown of the coalescence 
of two black holes. 

Until recently, numerical simulations were not stable and accurate enough for such 
estimates. This situation has changed dramatically in the past year, with several groups de- 
veloping codes that have allowed the evolution of binary black hole spacetimes from close to 
the ISCO through merger and ringdown. Most of these codes utilize techniques which allow 
the black holes, with their inhe rent singularities, to m ove successfully through the computa- 
tional domain. One approach (jPretoriusI 120051 . 120061 ) is to excise the singular region in the 
(physically inaccessible) interior of the black hole. In contrast, our metho d allows the singu- 



lar region to be approximately represented in the computational domain (IBaker et al.l 12005 



20061 : ICampanelli et al.l l2005bl : ICampanelli. Lousto. fc Zlochowerl 120061 ) . These new tech 



niques have led to dramatically more effective numerical simulations of binary black hole sys- 
tems, recently allowing us to determinate a ccurate waveforms for the final orbits and merger 
of equal mass systems (IBaker et al.l 120061 ) . which we generalize here for nonequal masses. 
While some preliminary n umerical calculations of recoil from mergers of non- spinning black 
holes have been reported (jCampanellil l2005at iHerrmann. Shoemaker. &: Lagunal 120061 ) , the 
initial separations have been too small and the resolutions too coarse for reliable and precise 
numbers. 

Here we report the first precise fully numerical estimates of the kick received from the 
merger of two nonspinning black holes. We choose a 1.5:1 ma ss ratio becaus e it is close to the 
analytically estimated optimal mass ratio for maximum kick (jFitchettlll983l ). but is also close 
enough to equal mass that resolution issues are not serious impediments to the numerical 
evolution. In § 2 we describe our numerical method and results, including convergence tests. 
We discuss the astrophysical implications of these results in § 3. 
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Numerical Simulations 



Reliable simulations of binary black hole mergers require the specification of constraint- 
satisfying initial field data followed by stable evolution of the Einstein equations. The evo- 
lution variables hold information about gravitational fields in the form of tensor fields rep- 
resenting the curvature of a vacuum spacetime. Einstein's equations, together with a crucial 
specification of gauge (i.e., coordinate) conditions, then govern the evolution of these fields 
in time. 



We utilize the "puncture" approach ((Brandt fc Brugmannlll997l ) to specify constraint- 
satisfying initial field configurations for black holes on an approximately circular inspiral 
trajectory within a few orbits of merger. Specifically, we consider three initial separations 
of coordinate distances d init = 4.1M , d init = 6.2M , and d init = 7.0M , where M is the 
total initial gravitational mass of the system, and we use units where Newton's gravitational 
constant G and the speed of light c are set to unity so that all quantities can be represented 
in terms of their mass-scaling. Further, we try several slight variations in the initial angular 
momentum for the largest separation. These data begin our simulations as the system 
approaches the ISCO, allowing about 1-3 orbits before the formation of a common horizon 
representing the final black hole. 



This data is evolved by our finite-differencing code, Hahndol (llmbiriba et al.ll2004l). o n 
an adaptive mesh refinement structure implemented via PARAMESH (iMacNeice et al. 
The punctures are allowed to move freely through the grid and are evolved according to a 
va riant of the BSSN f ormulation of Einstein's equations and gauge conditions as described 
in 



Baker et all ((20061); Ivan Meter et all ((2006). 



We interpret our numerical results by studying the gravitational radiation generated by 
the merger which carries gauge-invariant information away from the strong-field region. In 
our simulations, radiation is represented by a component of the space-time Weyl curvature 
tensor, ip±. In terms of t/u, the ti me rate of change of radiated momentum can be expressed 
as follows ((Newman fc Todll!980l ): 




We compute ip^ numerically, and extract it from the simulation data on a sphere of radius 
r = 50Mo, which we have found is sufficiently large enough to give precise results. 

Our main numerical results are given in Figured! which shows the net speed of the center 
of mass as a function of time for several runs. To compare runs at different initial coordinate 
separations we have shifted the time axis so that the kick speed peaks at t = 0. Note that 
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these curves are representative of a much larger body of simulation data. We have varied, 
for example, the finest grid spacing from hf — M /32 to hf — M /48 for the d-^ = 4.1M 
runs, and from hf — M /32 to hf = M /40 for the d init = 7.0M runs, and in either set 
of runs found the kick results of different resolutions agree to within 2%, thus confirming 
our numerical precision. (And w e have found that nume rical error in these runs converges 
away at an acceptable rate; see (IBaker et al.ll2005l . 120061 1 for previous convergence tests of 
our code.) We have also varied the extraction radius from r = 30M to r = 60M and, 
disallowing cases where the extraction radius significantly intersects refinement boundaries, 
we have found agreement across extraction radii to better than 2%. We have also varied the 
initial, orbital angular momentum in the d^t = 7.0M runs by as much as 4%, resulting in 
a 6% variation in the final kick values. Our full range of final kick values is 86-116 km s _1 . 
Excluding only the c? init = 4.1M case on the grounds that it starts with too little initial 
separation, our astrophysically relevant range of kick values is 86-97 km s -1 . 

The selection of our initial momenta and centers of mass was informed by previous work 



on quasicircular initial data based on minimization of an effec tive potential (jPfeifferl 12005 



Pfeiffer.Teukolskv & Cookll2000l ; Baker et al.ll2002l ; ICooklll994r ). and further refined by trial 
and error. The most astrophysically relevant simulations were presumed to be those with the 
smoothest, monotonic radiation frequencies and time derivatives thereof. By these criteria 
our most reliable run gave a final kick of 92 km s _1 . 

Also shown in Figure [T] is a post-Newtonian result for the kick, as integrated from 
a low- frequency cutoff. Because our simulation starts from a finite orbit, the resulting 
radiation has an effective low-frequency cutoff, which we have used in the post-Newtonian 
calculation for consisten cy. We integrate, specifically, a 2PN formula for the orbital frequency 
(IBlanchet et al. I Il995l ) and use the r esulting frequency as a function of time in a 2PN 
formula for the radiated momentum (IBlanchet. Qusailah. fc Willi 120051 ) . For most of the 
first orbit, the agreement of the = 7.0M simulation represented in Figure [T] with the 
post-Newtonian calculation is better than 1%. 

Recently, two groups have made refined analytic estimates of the kick from the merger 
of nonspi nning black holes, both with precis e answers but differing in magnitude by a factor 
of three. IBlanchet. Qusailah. fc Will J2005h pred ict a speed of 155±25 km s 1 for a 1.5:1 



mass ratio, whereas iDamour fc Gopakumarl (120061 ) predict 48 km s -1 . Our result is between 
these estimates but inconsistent with either one. We note, however, that Blanchet, Qusailah, 
& Will (2005) integrate only until the horizons overlap, which is close to the peak of our 
kick. At that point, our most reliable run (i.e., with the widest initial separation and least 
apparent eccentricity) gives a kick of « 150 km s _1 , which is encouragingly close to the 
analytic estimate. This points out the importance of the ~ 40% post-peak reduction in the 
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kick. This reduction can also be seen in Figure 5 of Damour & Gopakumar (2006). We 
attribute it to the evolution by more than tt radians of the phase of the emitted momentum 
during the merger, which thus partially opposes the vector kick that exists at the time of 
horizon overlap. It is also noteworthy that, integrating only up to ISCO, Blanchet, Qusailah, 
& Will (2005) obtain a kick of ~ 14 km s" 1 for a mass ratio of 1.5:1, and in our most reliable 
run we also obtain a kick of ~ 14 km s _1 at ISCO. 



An initial numerical estimate for 1.04:1 and 1.18:1 mass ratios was made by lHerrmann. Shoemaker, fc I 



(120061 ) (albeit at much lower resolution and starting much closer than our simulations), and 
they find kick speeds of 9 km s -1 and 33 km s -1 , respectively. One can compare results at 
different mass with reference to the lFitchettl (119831 ) fitting formula kick oc q 2 (q — l)/(l+g) 5 = 
(/j,/m) 2 dm/m, where q = m\jm<i > 1 is the mass ratio, fi/m = m 1 m2/(mi + m 2 ) 2 is the 
symmetric mass ratio and dm/m = {mi — m 2 )/(m 1 + m<i) is the fractional mass difference. 
Although originally derived in the context of a leading order post-Newtonian approximation, 
this f ormula is suggested to be a good approximation by perturbative Schwarzschild calcula- 
tions (IFitchett fc Detweilerlll984T). and it also closely agrees with more re cent post-Newtonian 



results (IBlanchet. Qusailah. Sz Willi 120051: iDamour Sz Gopakumar! 120061 ) . By Fitchett's for- 
mula, the llierrmann. Shoemaker, fc Laguna ( 2006 ) results would imply ~ 80 km s _1 for a 



mass ratio of 1.5:1, closer to our final value than the analytic estimates. 



3. Discussion and Conclusions 



Our range of 101 ± 15 km s _1 (and best estimate of 92 ± 6 km s _1 ) for the kick received 
in a 1.5:1 merger of nonspinning black holes has impor tant implications fo r the assembly of 
supermassive black holes in the early universe (see also iMerritt et al.ll2004l ). This is because 
as dark matter halos merge in the process of hierarchical structure assembly, their central 
black holes are also presumed to merge if the halo mass ra tio is not too extre me (otherwise 
dynamical friction on the halos is likely to be inefficient; see lTaffoni et al.ll2003l ). If black hole 
masses are linked to the mass of their host halos, it is therefore expected that in the early 
universe, black hole mergers are likely to involve comparable-mass objects. If the resulting 
kick exceeds the escape speed of the merged dark matter halo, the halo is left without a black 
hole. This could, therefore, have a significant impact on the number of black hole mergers 
in the early universe. In addition, if mergers between comparable-mass halos are common 
at redshifts z > 10, these ejections might reduce substantially the fraction of halos that host 
black holes. 



A full appraisal of the consequences of recoil will require detailed numerical studies of 
the effects of black hole spin and mass ratio, which will be our subject in future papers. 
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We can make initia 



gues ses by focusing on nonspinning black holes and by adopting as 



before the iFitchettl (119831 ) formula for the dependence of kick speed on mass ratio. Setting 
^kick(<7 — 1-5) = 92 km s _1 fixes the curve. 

To estimate th e escape speed from a dark matter halo, we follow the treatment of 
Merritt et al.l (120041 ) . The escape speed from a halo of virial mass M vir is V^ sc = 2cg(c)M vir / R vu 
where as before we set G = 1. Here c is the concentration parameter for the halo, g(c) = 
[ln(l + c) — c/(l + c)]~ , and we define R vir as the radius inside which t he average density 
is A r = 187T 2 + 82x — 39x 2 for a flat universe, where x = VIm{z) — 1 (see iBryan fc Norman 



1998 



equation 6). The critical dens ity for a flat un iverse is p cr it = 3ffo/(87r)[pM(l + z ) 3 + 



Ov] 
tt A = 

as c 



10 -29 g cm -3^ + ^3 



e.g 



0.73. ^From equation (18) of iBullock et al.l (120011 ). the concentration parameter scales 



Peacock! Il999l. sec tion 3.2). We use Qm = 0.27 and 



9[M vir /(2xlO 13 M )]- ai3 (l + ^) 



In Figure[2]we plot the minimum halo mass as a function of redshift such that V^ sc > Wkick 
for the listed mass ratios of 1.5:1, 3:1, 5:1, and 10:1, as projected from our results by the 
Fitchett scaling. Note that the low power of M vir in the escape speed means that the 
minimum halo mass depends sensitively on the kick speed. For example, if the kick speed at 
a 1.5:1 mass ratio were 200 km s _1 instead of 92 km s , the threshold halo mass at z — 10 
for this mass ratio would jump from 4 x 10 8 M Q to 3 x 1O 9 M . This would in turn reduce 
the number density of halos massive enough to retain their black holes by a factor ~ 20 (see 
Mo fc White! 120021 for a pedagogical discussion of how to estimate halo number densities). 
This points out the importance of estimating kick speeds accurately. 



As discussed by iMerritt et al.l (120041 ). kicks are also important in the current-day uni- 
verse for low-mass concentrations of stars. Their Figure 2 is a useful summary of escape 
speeds from the centers of galaxies and globular clusters. From their figure, we see that 
comparable-mass mergers with kick speeds ~ 100 km s _1 will cause ejection from globulars 
and some dwarf ga laxies, but that more massive galaxies will retain the remnant. Using 
the IFitchettl (119831 ) scaling with mass ratio, we find that black holes that are > 10 times 
more massive than their companions are not ejected from globular clusters with escape 
speeds of a few tens of km s" 1 . Coincidentally, this is approximately the same mass ratio 
required to protect a massive black hole fr om cluster ejection from three-body interactions 
(jGiiltekin. Miller, fc Hamilton! 120041 . 120061 ) . Therefore, if an intermediate-mass black hole 
of > 100 — 200 M Q is formed in a cluster, it can stay and potentially grow through future 
interactions. Stellar-mass binary black holes (M < 50 M m ) will be ejected from such clusters 



prior to merger by Newtonia n three-body interactions (IKulkarni. Hut, fc McMillan! 11993 



Sigurdsson &: Hernquistlll993l ). hence recoil is not so important for low- mass black holes in 
this context. However, mergers of black hole binaries in low-density galactic disks could 
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produce a population of high-speed coalesced black holes. 

In conclusion, we have presented a reliable and precise, fully numerical estimate of the 
gravitational recoil produced by the merger of two unequal mass nonspinning black holes. 
Our best estimate of 92±6 km s" 1 for a 1.5:1 mass ratio is intermediate between the recent 



analytic estima tes oflBlanchet. Qusailah. fc Willi (120051 ). who suggest 155±25 km s 1 for this 



mass ratio, and lDamour fc Gopakumarl (120061 ). whose formulae would imply 48 km s 1 . Our 



results are thus an important step in accurately evaluating the astrophysical consequences 
of gravitational radiation recoil in dense stellar clusters and the early universe. 

We gratefully acknowledge supercomputing support from the Commodity Cluster Com- 
puting Project (NASA-GSFC), the NASA Center for Computational Sciences (NASA-GSFC) 
and Project Columbia (NASA Advanced Supercomputing Division, NASA- Ames). This work 
was supported in part by NASA grants ATP02-0043-0056 and O5-BEFS-05-0044. JvM, MK, 
and MCM were also supported in part by the Research Associateship Programs Office of 
the National Research Council and Oak Ridge Associated Universities. Comments by Scott 
Hughes and Cliff Will on a earlier draft helped clarify this manuscript. 



REFERENCES 

Baker, J., Centrella, J., Choi, D.-L, Koppitz, M., & van Meter, J. 2006, Phys. Rev. Lett., 
96, 111102 

Baker, J., Centrella, J., Choi, D.-L, Koppitz, M., & van Meter, J. 2006, Phys. Rev. D, 73, 
104002 

Baker, J., Campanelli, M., Lousto, C. O. & Takahashi, R. 2002, Phys. Rev. D, 65, 124012 
Bekenstein, J. D. 1973, ApJ, 183, 657 

Blanchet, L., Damour, T., Iyer, B. R., Will, C. M., & Wiseman, A. G. 1995, Phys. Rev. 
Lett., 74, 3515 

Blanchet, L., Qusailah, M. S. S., & Will, C. M. 2005, ApJ, 635, 508 
Boylan-Kolchin, M., Ma, C.-P, & Quataert, E. 2004, ApJ, 613, L37 
Brandt, S. & Briigmann, B. 1997, Phys. Rev. Lett., 78, 3606 
Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80 



- 8- 



Bullock, J. S., et al. 2001, MNRAS, 321, 559 
Campanelli, M. 2005, Class. Quant. Grav., 22, S387 

Campanelli, M., Lousto, C. O., Marronetti, P., & Zlochower, Y. 2006, Phys. Rev. Lett., 96, 
111101 

Campanelli, M., Lousto, C. O., & Zlochower, Y. 2006, Phys. Rev. D, 73, 061501 

Cook, G. B. 1994, Phys. Rev. D, 50, 5025 

Damour, T., & Gopakumar, A. 2006, Phys. Rev. D, 73, 124006 

Favata, M., Hughes, S. A., & Holz, D. E. 2004, ApJ, 607, L5 

Fitchett, M. J. 1983, MNRAS, 203, 1049 

Fitchett, M. J., & Detweiler, S. 1984, MNRAS, 211, 933 

Giiltekin, K., Miller, M. C, & Hamilton, D. P. 2004, ApJ, 616, 221 

Gultekin, K., Miller, M. C, & Hamilton, D. P. 2006, ApJ, 640, 156 

Haiman, Z. 2004, ApJ, 613, 36 

Herrmann, F., Shoemaker, D., & Laguna, P. 2006, |gr-qc/0 601026 

Imbiriba, B., et al. 2004, Phys. Rev., D70, 124025 

Kulkarni, S. R., Hut, P., & McMillan, S. L. W. 1993, Nature, 364, 421 

Libeskind, N. I., Cole, S., Frenk, C. S. & Helly, J. C. 2005, MNRAS, 368, 1381 

MacNeice, P., Olson, K. M., Mobarry, C, deFainchtein, R. & Packer 2000, Computer Physics 
Comm., 126, 330 

Madau, P., & Quataert, E. 2004, ApJ, 606, L17 

Merritt, D., Milosavljevic, M., Favata, M., Hughes, S. A., & Holz, D. E. 2004, ApJ, 607, L9 



Micic, M., Abel, T., & Sigurdsson, S. 2005, MNRAS, submitted (jastro-ph/0512123) 
Miller, M. C, & Colbert, E. J. M. 2004, IJMPD, 13, 1 
Miller, M. C, & Hamilton, D. P. 2002a, MNRAS, 330, 232 
Miller, M. C, & Hamilton, D. P. 2002b, ApJ, 576, 894 



- 9- 

Mo, H. J., & White, S. D. M. 2002, MNRAS, 336, 112 

Mouri, H., & Taniguchi, Y. 2002a, ApJ, 566, L17 

Mouri, H., & Taniguchi, Y. 2002b, ApJ, 580, 844 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Newman, E. T. & Tod, K. P. 1980, General Relativity and Gravitation: One hundred years 
after the birth of Albert Einstein, vol. 2, Ed. by A. Held, Plenum Press N. Y., (1980), 
pp. 1-36 

O'Leary, R. M., Rasio, F. A., Fregeau, J. M., Ivanova, N., & O'Shaughnessy, R. 2006, ApJ, 
637, 937 

Peacock, J. A. 1999, Cosmological Physics (Cambridge: Cambridge University Press) 
Peres, A. 1962, Phys. Rev., 128, 2471 
Pfeiffer, H. P. 2005 ( |gr-qc/05 10016ft 

Pfeiffer, H. P., Teukolsky, S. A., Cook, G. B. 2000, Phys. Rev. D, 62, 104018 
Pretorius, F. 2005 Phys. Rev. Lett., 95, 121101 
Pretorius, F. 2006 ( gr-qc/0602115| 

Redmount, I. H., & Rees, M. J. 1989, Commun. Astrophys., 14, 165 

Sigurdsson, S., & Hernquist, L. 1993, Nature, 364, 423 

Taffoni, G., Mayer, L., Colpi, M., & Governato, F. 2003, ApJ, 341, 434 

Taniguchi, Y., Shioya, Y., Tsuru, T. G., & Ikeuchi, S. 2000, PASJ, 52, 533 

van Meter, J. R., Baker, J. G., Koppitz, M. & Choi, D.-I. 2006, Phys. Rev. D, 73, 124011 

Volonteri, M., & Perna, R. 2005, MNRAS, 358, 913 

Wiseman, A. G. 1992, PRD, 46, 1517 

Yoo, J., & Miralda-Escude, J. 2004, ApJ, 614, L25 



This preprint was prepared with the AAS IATgX macros v5.2. 



-10- 



180 




time (M) 

Fig. 1. — The magnitude of the radiated momentum, as a function of time, from three 
different simulations. For initial coordinate separations of d m \ t = 4.1M , 6.2M , and 7.0M , 
the final values of the kicks are respectively 113 km/s, 97 km/s, and 92 km/s. Also shown 
is the 2nd order post-Newtonian radiated momentum, which was computed from a low 
frequency cutoff commensurate with that of the d^it = 7.0M simulation (see text for details). 
The excellent agreement of the post-Newtonian kick with that of the dj n j t = 7.0M simulation 
over most of the first orbit, together with the agreement to within 6% of the final kick from 
the cfjnit = 7.0M simulation with that of the di^t = 6.2M simulation, lends support to the 
accuracy of these results. 
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Fig. 2. — Minimum mass of a dark matter halo at a given redshift required to retain the 
product of the merger of two nonspinning black holes with a mass r atio indic a ted on the curve. 
Details of the computation are in the text. Note that we use the iFitchettl (I1983I ) analytical 
estimate of the mass ratio dependence, assuming our numerical result of 92 km s _1 for the 
mass ratio of 1.5:1. Then a mass ratio of 3:1 gives 140 km s _1 , a mass ratio of 5:1 gives 
103 km s _1 , and a mass ratio of 10:1 gives 45 km s _1 . This figure indicates that early halos 
might lose merger remnants because of the kick, but in the current universe only globular 
clusters or the smallest dwarf galaxies could have black holes ejected because of gravitational 
radiation recoil. 



